RIKEN-QHP-45, RIKEN-MP-55 



Relativistic hydrodynamics from projection operator method 

Yuki Minami and Yoshimasa Hidaka 
Theoretical Research Division, Nishma Center, RIKEN, Wako 351-0198, Japan 

(Dated: October 5, 2012) 

We study relativistic hydrodynamics in the linear regime, based on Mori's projection operator 
method. In relativistic hydrodynamics, it is considered that ambiguity about the fluid velocity occurs 
from a choice of a local rest frame: the Landau and Eckart frames. We find that the difference of 
the frames is not the choice of the local rest frame, but rather that of dynamic variables in the linear 
regime. We derive hydrodynamic equations in the both frames by the projection operator method. 
We show that natural derivation gives the linearized Landau equation. Also, we find that, even for 
the Eckart frame, the slow dynamics is actually described by the dynamic variables for the Landau 
frame. 
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I. INTRODUCTION 



Relativistic hydrodynamics has been widely apphed for studying relativistic nonequihbrium phenomena. For ex- 
amples, it describes hadron spectra and elliptic flow in the heavy ion physics [1, 2], and jets like gamma-ray bursts in 
the astrophysics [3, 4]. The hydrodynamic equations applied to these systems are mainly those for perfect fluids. One 
of reasons is that the dissipative effects in relativistic hydrodynamics are not fully understood, e.g, some pathological 
problems arise from the dissipative effects: the acausal propagation and the instability of the equilibrium state [5]. 
Although many hydrodynamic equations have been proposed to resolve these problems [6-11], it is not still obvious 
which equation describes the correct behavior of the relativistic dissipative fluid. Namely, even the basic equation has 
not been established in relativistic hydrodynamics. 

The relativistic hydrodynamic equations are generally given as the following conservation laws: 

df.r = 0, (1) 

df.T'"' = 0. (2) 

Here, j'^ is the particle current and T'"^ is the energy-momentum tensor. They are decomposed to 

/' = nu^ + z^'', (3) 
T^"' = hu^'u" - Pg^''' + q^'u" + q^u^' + r^'', (4) 

where n is the particle density, h = e + P the enthalpy density, P the pressure, e the energy density, the fluid-four 
velocity. The dissipative terms, i^^, q^, and r^"^, denote the particle and energy diffusions, and the viscous stress 
tensor, respectively. The explicit expressions of these terms are not unique but depending on considered equations. 
This ambiguity comes from a choice of local rest frames of the fluid. 

To see this ambiguity, let us classify the hydrodynamic equations into two groups: the Eckart and Landau frames [12, 
13]. In the Eckart frame, the local-rest frame is defined as that of the particle current, i.e., the fluid velocity is 
proportional to the particle current: 

u^cxf. (5) 

In this frame, the particle diffusion is absent, v'^ = 0. On the other hand, in the Landau frame, the fluid velocity is 
proportional to the energy current: 

< cx <T,^ (6) 

In contrast to the Eckart frame, the energy diffusion is absent, = 0. We note that nonrelativistic hydrodynamics do 
not have these ambiguity. In the nonrelativistic limit, the energy current is identical to the particle current because 
the mass energy dominates the energy of fluids. Actually, the Navier-Stokes equation does not have such ambiguity 
and the frames. It is considered that this difference between the frames is just by the references frames and apparent. 
However, several differences, which are not just apparent, actually exist. For example, the Eckart frame has the 
instability of the equilibrium state, but the Landau frame does not. 

To discuss the difference of the Landau and Eckart frames, we consider fluctuations from the global equilibrium 
state, namely, the linear nonequihbrium regime. The merit of this fluctuating state is that we can observe the state 
at the same rest frame for the energy and particle currents. We note that, at the equilibrium state, the particle and 
energy currents rest: = = (1,0). Then, in the fluctuating state, we also have the same reference frame for 
the Landau and Eckart frames because the considered state is just perturbed from the equilibrium one; moreover, we 
need not to bother about what are local equilibrium and local rest for the relativistic system. 

To see relativistic hydrodynamics in the linear regime, let us consider the Landau and Eckart equations as examples. 
For the Landau equation, the dissipative terms read 

^^ = a(^^)V(/3m), (7) 

= 0, (8) 
2 
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-CA^^^id^-u), (9) 



where A, rj and ( are the thermal conductivity, the share and bulk viscosities, respectively. A^'' = g'"^ — u'^u'^ is a 
projection and = A.^^'^d^, is the space-like derivative. 
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Now, wG linearize the Landau equation about fluctuations from the equihbrium state. Let us write n(x) = no+5n(x), 
e{x) = Co + 6e{x), P{x) = Pq + SP{.x), {l3^i)(x) = (/3//)o + S{l3^i){x), and u^'{x) = + 6u^[x). Here, the symbols 
with the prefix 5, denote the fluctuations. The equilibrium values are denoted by the sufiix 0. Hereafter, we employ 
variables with the suffix and the prefix as the equilibrium values and fluctuations, respectively. For simplicity, let us 
choose the rest frame as the reference frame: Uq = (1, 0). Then, by the relation in the linear regime, UgSuf^ = 0, the 



fluid-velocity fluctuation is written as 

Suf"^ i0,6vi^). (10) 
In consequence, the Landau equation is linearized as 

dodn=-noV-6vi^ + x(^^^ V^S{Pfi), (11) 

doSe ^ -hoV ■ 5vl, (12) 

doihoSvi^) = -V{SP) +(^C+ V(V ■ Svi^) + ijV^Svi^. (13) 

Next, let us move on to the Eckart equation, the dissipative terms of the Eckart equation are 

1^" = 0, (14) 

qf" = \{dlT -TDu^'), (15) 

T^"' = 7] df^u" + dluf" - ^A'^'id^-u) + (A'^'id^-u), (16) 

o 

where D = u^d^ is the time-like derivative. The linearized equations are 

doSn = — noV • Sve, (17) 

do5e = -hoV ■ Sve + X{V^6T + TodoV ■ Sv^), (18) 

do{ho6vE) = -V{SP) + r^VHvE + f C + ^»?Vc^ " ^^^) + ^{^do6T + T^dlSv^ . (19) 



We note that the linearized Landau and Eckart equations have different forms even in the same rest frame, as we 
mentioned. 

To investigate relativistic hydrodynamics in the linear regime, we here use Mori's projection operator method [15]. 
Mori's projection operator method is a powerful tool for extracting slow dynamics. This method is widely applied 
and succeed in condensed matter physics [16-18]. Actually, various slow dynamics, e.g. the Navier-Stokes, Langevin, 
Boltzmann equations, and equations for Nambu-Goldstone bosons, are derived [17, 19, 20]. The merit of the projection 
operator method is that we can derive slow dynamics only by choosing slow variables and commutation relations of 
those without microscopic details. We note that dynamics on macroscopic scale can be described by much fewer 
degrees of freedom than those on microscopic scale. Such degrees of freedom are called the slow variables (or gross 
variables). The slow variables are degrees of freedom that label a macroscopic state and describe long-time behavior. 

The view point of the projection operator method tells us that a choice of the slow variables is a key for hydrody- 
namics. From the linearized Landau and Eckart equations, (11)-(13) and (17)-(19), we see that the dynamic variables 
are given as the energy and particle densities, and fluid- velocity fluctuations [33]: 

{5e,5n,5vl,,}. (20) 

The fluid velocities are different depending on the frames: 

<54 = n-^j\ (21) 
K = V^T""- (22) 

Namely, the difference of the frames is the choices of the slow variables. 

The important point about the choices is that f is not essentially slow because it is not a conserved charge. In 
contrast, T^^ is a conserved charge and slow because the energy current is equivalent to the momentum for the 
relativistic system: T*^' = T'°. Actually, we shall find that a slow part in originates from T^'. We will provide 
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detail on this point in Sec. II C and V. Then, in this paper, we will apply the projection operator method for the 
choices, 

{5e,5p\6n}, (23) 
{Se,Sp\Sn,Sf}, (24) 

where Sp^ = and Sj'^ = f . We will show that the first set of slow variables, Eq. (23), gives the linearized Landau 
equations (11)-(13). Furthermore, we will derive the equations for the second set, Eq. (24), which include the Landau 
and Eckart frames. Then, we will derive the linearized Eckart equations by eliminating Sp^ from the equations for 
{5e, (5n, (5j*}. We note that, to correctly treat the slow part of 5j'^ coming from we have to choose the both 
of them as the slow variables. 

After that, we will discuss properties of the derived equations. We shall show that the equations for {(5e, (5p', (5n} 
and {(5e, (5p*, Jj*} have the same slow modes. Namely, Sj"^ contains the irrelevant part for the slow dynamics 
of relativistic fluids. Moreover, we will illustrate that, even for the Eckart equation, the slow dynamics is actually 
described by the Landau's variables, {Se^Sp^ ,5n\. 

Our study is first attempt for the derivation based on the projection operator method. The earlier studies [7-10] 
assume that the relativistic Boltzmann equation as an underlying microscopic theory. In contrast, we stress that our 
derivation is independent of microscopic details. 

This paper is organized as follows. In Sec. II, we briefly review Mori's projection operator method for readers 
being unfamiliar with it. Also, we explain conserved charges as the slow variables and discuss those for relativistic 
hydrodynamics. In Sec. Ill, we first determine equal-time correlations of the slow variables to determine the property 
of the equilibrium state. We note that we here consider the fluctuations from the equilibrium state. Our determination 
is based on thermodynamics and Lorcntz symmetry on microscopic scale. Then, it is independent of the microscopic 
detail. In Sec. IV, we derive slow dynamics by the projection operator method for the sets, Eqs. (23) and (24). In 
Sec. V, we discuss the results obtained in Sec. IV. In particular, we study slow modes of the equations for {Se, Sp^, 6n} 
and {5e,Sp^,6n,Sj^}. We explicitly show that these equations have the same modes. Furthermore, we consider the 
Onsager reciprocal relation in the Eckart equation to illustrate that the slow part of Sj coming from {^e, Sp^, Sn}. 



II. MORI'S PROJECTION OPERATOR METHOD 



In this section, we provide Mori's projection operator method [15, 17, 18, 21, 22]. We can formally extract slow 
dynamics from microscopic Hamiltonian dynamics by this method. On the microscopic scale, an operator at time t, 
0{t) = 6*^*0(0)6"*^*, evolves by the Heisenberg equation, 

doO{t) = i[H, 0{t)] = iCO{t), (25) 

where C is the Liouville operator. In the following, we decompose this time evolution equation into slow and fast 
ones. 



A. Projection operator 

First, we introduce basic ingredients for the projection operator. Let us consider a many-body system at finite 
temperature. As an equilibrium distribution, we assume the grand-canonical one. Then, the density matrix is given 

as 

g-/3(ff-MAf) 

= tre-/3(^^-M^)' ^^^^ 

where H is the Hamiltonian, N is the number operator, /x the chemical potential, and the inverse temperature /3 = 1 /T. 
With the density matrix, thermal average of 0{t) is defined as 

(0(t))eq = tr PeqO(t) = (O(0))eq. (27) 

Also, we define an inner product of A and B as 



P Jo P Jo 



(28) 
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Moreover, we introduce a set of slowly-varying operators (slow variables), {An(t,x)} = {Ai,A2, A„}. If we can 
separate the time scale into long- and short-time ones, such operators exist and describe the slow dynamics. Let us 
consider the slow operators at the initial time t ~ 0, {An{0,x)}. In general, they arc not orthogonal to each other. 
We introduce metric to consider the orthogonal basis: 

gnm{x - y) = (A„(0, x),Arn{0, y)). (29) 

The orthogonal operators, represented with an upper index, is defined as 

A"{t,x) = / d'yg'"'\x - y)A^{t,y), (30) 



where g^"^{x — y) is the inverse of gnm{x — y). These quantities are orthogonal to those with lower indices, 

(A„(0, x), A™(0, y)) = 5.:^5{x - y), (31) 
E / d'^ygn^nix - y)g"''{y - z) = S,^5(x - y). (32) 

It is useful to see the relation between the metric and the effective potential V^({^„}), which is given as 

^({^4) = ^ E / dxdyA,{x)g"^^\x - y)AM + (33) 

n,m 

where ... are higher order terms of A„ and can be neglected in the linear regime. We note that we are interested in 
fluctuations as A„. 

We have prepared the basic ingredients. Let us introduce the projection operator acting on any operators B{t, x) 

as 

VB{t,x) = E / rf'y^n(0,y)(B(t,y), A"(0,a;)). (34) 

n •' 

The projection operator extracts the slowly varying part of B, which is determined only by the slow variables {A„}. 
We also define the orthogonal projector as Q = 1 — for later use. 

B. Generalized Langevin equation 

In this subsection, we derive so-called the generalized Langevin equation. This equation is given by decomposing 
the Heisenberg equation into slow and fast parts. For the decomposition, we use the following operator identity: 



Jo 



(35) 



which is valid for arbitrary C and V [16]. 

Let us derive this identity. First, we consider the following decomposition: 



age'^* = e'^HC = e'^'ViC + e'^* QiC. (36) 
Next, we consider the Laplace transform of exp{iCt), 



' 1 

dte-'^e'^* = -. (37) 







Then, we decompose Eq. (37) into 



^ ^ :iz-QiC)- ^ 



z ~ iC z — iC z — QiC 

^ -,{z-iC + ViC) ^ 



z — iC z — QiC 
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Performing the inverse Laplace transform, we find the identity, 

e''^' = e^''^^ + f dse'^'^*-'^PiCe^'^'. (39) 
Jo 

Substituting Eq. (39) into the second term of Eq. (36), we obtain the operator identity, Eq. (35). 

Muhiplying Eq. (35) by an initial value of the slow operator, we obtain the decomposed equation of motion for 
A„{t) = e*^*A„(0): 



doAr,{t,x) = J (fyinn'^'ix - y)A„,{t,y) - dscPv<Pn"\t - s,x - y)A„,{s,v) + Rn{t,x), 
without any approximations. Here, we introduced the following functions and operator: 



(40) 



iQrTix -y) = (z£A.(0, x),A''\Q, y)) = -^*(K(0, x), A'"t(o, y)])eq, (41) 

^,r{t-s,x-y) = -0{t-s){t£R^{t,x),A"'{s,y)), (42) 
i?„(i, x) EE e''^^QiCAn{0, x). (43) 

Equation (40) is the generalized Langevin equation and has the following properties: 

1. Equation (40) is the operator identity. 

2. The first and second terms in the right-hand side represent the slow motions. 

3. The first term corresponds to a time-reversible change. 

4. The second term corresponds to a time- irreversible change. Also, this term depends on a past time value, Am{s), 
for s < t. <?„™(t — s,x — y) is called the memory function. 

5. The last term is the noise term corresponding to the fast motion. For hydrodynamics, this term is usually 
neglected. On the other hand, for the Langevin dynamics, we treat this term as a random noise. 

It is useful to rewrite Eq. (40) to the equation in momentum space: 

doAn{t, k) = inn"\k)A,n{t, fe) - / ds<?„™(i - S, fc)A«(s, fe) + Rn{t, k) . (44) 

Jo 

For time component, we perform the Laplace transform: 

yl„(z,fc) = j dt j d^xe- A„{t,x). (45) 

Then, Eq. (44) becomes 

zA{z, k) = iQ{k)A{z, k) - <P(z, fc)A(z, k) 4- R{z, k) + A{t = 0, k) (46) 
in the Laplace-momentum space. Here, A[t = 0, k) is the initial value and we used matrix notation. 

C. Conserved charges as slow variables 

Here, we provide why conserved charges are slow and discuss the dynamic variables for the Landau and Eckart 
frames. The key is that conserved charge densities generally satisfy conservation laws: 

doj' = ~d,j\ (47) 

where is a conserved charged density, and is its current. In the momentum space, the conservation law becomes 

do.f = tk,f. (48) 

We note that the time change rate of j*^ is proportional to the wavenumber, so that the low-wavenumber components 
turn out to be slow. Therefore, the change of the conserved charge densities are necessarily slow in the low-wavenumber 
region, i.e., on macroscopic scale. 
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Now, let us consider the case of relativistic hydrodynamics. We here have the three conservation laws in Eqs. (1) 
and (2). From those, we obtain the three conserved charges, the particle number, the energy and the momentum: 

dof = thf, (49) 
aoT™ = thT°\ (50) 
doT'° = ikjT'^ . (51) 

We note that the above quantities are the slow variables for the Landau equation. The important point is that the 
particle current is not conserved: 

dof^ikjU^'. (52) 

Thus, the time change rate is not proportional to the wavcnumber. Namely, the particle current is not essentially 
slow although it is proportional to the fluid velocity for the Eckart frame. 

Nevertheless, we note that the particle current has a slow part coming from the conserved charges, because is 
not orthogonal to {j", T^", T"*}. In other word, the projection of j' on those does not vanish, 

Vf + 0, (53) 

and gives the slow part. From this slow part, we can derive the linearized Eckart equation as we will show in Sec. IV. 
Here, we stress that the slow dynamics is essentially determined by the conserved charges, T"", T°'}, even for the 
Eckart equation. 



III. METRIC AND THERMODYNAMIC QUANTITIES 

In this section, we discuss relations between the metric (7„,„ and thermodynamic quantities [16, 17, 23]. As discussed 
in Sec. I, we employ the fluctuations of conserved charges as slow variables, i.e., A„ = {(5e, Jj*} with be = 

— eo, <5p* = bn = j° — no, and bf = with cq = (T°°)oq and uq = (n)eq. We assume that the density 
matrix at thermal equilibrium is invariant under time reversal transformation, i.e., l'peqT~^ = Peqi where T is the 
time reversal operator. The slow variables transform under T as 

TSe{t, x)T^^ = 5e{-t, x), r5n{t, x)T^^ = 5n{-t, x), (54) 
rSp'{t,x)T-^ = -5p\-t,x), TSfit,x)r-^ = -5j\-t,x). (55) 

5e and 5n (Sp^ and Sj^) are even (odd) operators, so that dn and Se does not mix Sp^ and Sj^, i.e., gep{k) = 9ej{k) = 
97ip{k) = gnj{k) = 0. 

Since we are interested in the low energy behavior of slow variables, we apply the derivative expansion. The metric 
is expanded as a power series of 

9nm{k) - gnm + git^k^ + g^nl,^,k'k^ + " ' " , (56) 

where we assumed gnm{k) is analytic at fc = 0; in other words, there are no long range correlations. The only 
leading terms, gnm, contribute to the linearized liydrodynamic equations at first order, so that we will not consider 
contributions from the higher order terms. 

A. Qnn^ Qee Jllici Qen 

First, we focus on (;„„, g^e and gen, which are given as 

gee = (fx{Se{x),Se{0)), gnn = x{6n{x) , Sn{0)) , = gne = (fx{Se{x),Sn{0)). (57) 



By definition, these metric are related the susceptibilities. In fact, gnn can be expressed using the isothermal com- 
pressibility kT as 

^ (((5^)2),q = n^Tofcr, (58) 



ToV 
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where 6N = N ~ {N)cq and V is the volume, gee and gen are related to the specific heat at constant volume 
cy = /3o(ffee ~ 9ne/9nn) = /^o/ff'^'^- The invcrse matrices for e and n can be written by the entropy density as a 
functional of e and n: 

9"' = -% = 5"^ = 5^" = -|5f, (59) 

ae^ an^ oeon 

because the effective potential as a functional of 5n and (5e is given as the entropy density s{5e, Sn), for grand canonical 
ensembles [16]. Then, the slow variables with upper index are written as 



If we notice the thermodynamic relations 



ds\ „ / ds 



we see that 



A^--[^)se-(^)sn^^Sf^, (63) 



de J \dn 



Namely, —d/3 and S{/3iJ,) are orthonormal to (5e and 6n: 



{5e{t, r), J(/3/.)(t, r')) = (<5n(i, r), <5/3(t, r')) = 0, (65) 
-((5e(t, r), <5^(t, r')) - {5n{t, r),Si(3t,)it, r')) - 5(r - r'). (66) 



B. gp^pj , gp^Jj , and g 



Next, we consider gpipj , (/pijj and ffji^j • We show that two of them, gpipj and gpijj , are expressed as the enthalpy 
and the number density: 



gp,p, = J d^x{T'^H0,x),T°H0,O)) = (67) 
gp,p = I d3a;(r«(0,a;),/(0,0)) = <5'^Tono. (68) 



The same relations are obtained in Ref. [23]. These relations are derived from Lorentz symmetry underlying theory. 
For an arbitrary Hermitian operator O, the following identity is satisfied: 



d''x(T"'(0, a;), O) = {P\0) = i{[H, K\ O) = -iT{[K\ 0])^q, (69) 

where is the boost operator, [H, K"^] = ~iP^ and the following Kubo's identity was employed: 

([H,A],B) = -r([A,i?t]),q. (70) 

Since T^^{x) and j^{x) are Lorentz tensor and vector, respectively, these transform under Lorentz transformation as 

[i'"^, T^P{x)] = i{x^'^'' - x''d^')T^P{x) - iiif^T^Pix) - ^f^T^'P{x) + r^^PT^" - ^fPT^^'{x)), (71) 
[L^^\3P{x)] = z(x'^9'' - x''9'^)/(x) - *(rr/(x) - r?''''j'^(a;)), (72) 
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where L'^'^ is the charge of Lorentz symmetry, and 77^"^ = diag(l, —1, —1, —1) is the (inverse) Minkowski metric. For 
the Lorentz boost = L'°, they obey 

[K\ T°^{x)] = ~i{x"d' ~ x'd°)T°^{x) + iT'^ ~ i?fT°"{x), (73) 

[K\f{x)] ^ -i(x°9* - x'd°)f{x) - i7j'^?i{x). (74) 

Therefore, the thermal averages for these commutators satisfy 

\ T"J»])eq = i{T''{x) - 77'JTO"(x))eq - iS^'ho, (75) 

{[K\f^{x)]U = -7(?/^^n(x))eq = tS^'no. (76) 

Inserting Eqs. (75) and (76) into Eq. (69), we arrive at Eqs. (67) and (68). These identities enable us to relate 
two-point functions to one-point functions. 



IV. APPLICATION OF MORI'S PROJECTION OPERATOR METHOD TO RELATIVISTIC 

HYDRODYNAMICS 



In this section, we apply Mori's projection operator method to relativistic hydrodynamic systems, and derive 
equations of motion for {Se, dp\ 5n} and {6e, 5p^, 6n, 5j'}. We first show that the set, {Se, Sp'', Sn}, gives the linearized 
Landau equation. For the Eckart equations, we introduce the current of the conserved charge, Sj'', which is proportional 
to the fluid velocity in the Eckart frame, in addition to 6p'' and Sn. We employ the derivative expansion and keep the 
spatial and time derivative to the second order, i.e., i9o, V, V^, i9o V and Oq. We will drop the noise term i?„(t, x) in 
the equation of motion. This term is irrelevant in the time evolution of the expectation value. If one is interested in 
stochastic hydrodynamics, one may keep the noise term [24]. 



A. Linearized Landau equation 



First, we derive the linearized Landau equation. For this purpose, we choose Se, Sp and 6n as slow variables. Since 
Sp^ is chosen as a slow variable, the equation for doSe does not contain dissipative terms, 

dQSe^-V-6p = -hoV -Svl, (77) 

where wc defined the fluid velocity Svl = Sp/h^. This equation is nothing but energy conservation law, Eq. (12). 
This can be confirmed by the the following calculation: 



if2P{k)= I d^xe-'''-''{iC6e{x),6p>{0))gP'p\k) 

(78) 



d^xe-'''-''{~\/^p^{x),Sp'{0))gP'P\k) 



-ik'gp.p,{k)gP'P (fc) 
-ik\ 



Therefore, the reversible term becomes —V • 6p. The memory function vanishes because iC6e turns out to be —ik ■ 6p, 
and then QiCSe = [see Eqs. (42) and (43)]. 

Let us move onto the equation for dgSn. For the reversible part, the only ifinp' survives in f2 from time reversal 
symmetry, which is 

in,,p,{k)^ I d^xe-''' ''{iCSn{x),Sp'{0)) 

= -ik^gj.p. + 0{k^) 
= -ik'Tonn + 0{k^). 
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For the memory function, we keep it up to of order k^. <?„e vanishes as the same as previous due to QiCSe = 0. 'Pnp 
is of order k'^ from tensor structure, which can be neglected. Therefore, we may only consider Since we are 

interested in slow dynamics, we also expand the memory function in terms of z in addition to k: 

<?„„(z, fc) = ^ dtj d^xe-''e-''' ''{e''^^QiC6n{0, x), iCSn{0, 0)) 

= / dt f d^xe-'''e-'''''{e''^^QSf{0,x),dj^{0,O)) 



fOO n (80) 

~fcM d^x{6f it,x)-'^Sp'{t,x), Sf (0,0)- '^6p\0,O)) 

Jo J ho ho 

where ~ denotes the approximation of order k^ and z". The approximation of order corresponds to the Markov 
approximation, i.e., in the coordinate space, <?„„(t, a;) ~ — X'V^ S{t)6^^^ (x) . We defined the thermal conductivity A as 



\noT(_ 



/ dt / d^x{5f{t, x) - Sf{0, 0) - ^dp\0, 0)), (81) 

Jo J ho tio 



and we used 



Q5f{t,x) 5f{t,x) - ^Sf{t,x) (82) 
ho 

in the leading order of the derivative expansion. We note that the second term is important to remove the contribution 
of zero mode from Sj^. As a result, we arrive at the equation for doSn as 

do6n = -noV • Svl + ~XV^6{f3ti). (83) 

This equation coincides with Eq. (11). 

Similarly, for 6p^ , the reversible term if2pe = —ik'^Toho, ifipn = —ik^ToUo, so that 

-iHpeFe - lOpnFn = -ik^-TohoSji + ToJiQ^ (/^aO) = -ik'SP, (84) 
where we used Gibbs-Duhem's relation, 

5P^^ST + Tono5{P^x). (85) 
Jo 

For the dissipativc terms, 4>j,e vanishes, and ^ k^ from tensor structure can be neglected as the same as the 
previous. Therefore, only fppp survives in the leading order, which is evaluated as 

^>p.pfc(z,fc) ~ Pfc' I dt j dx{T\t,x),T^\Q,x)) 

^ (86) 
= To(c+^vYk^+ToVk^5'\ 

where we used the same approximation in Eq. (80). The shear and bulk viscosities arc defined by Kubo formula as 



V = P^j^ dtj d'xiT'^{t,x),T'\0,O)), (87) 
C-|r; = /?o^°°rfty d^x{T^\t,x),T^^0,O)). 



Noting that = —/3oSv)^ = —{l3o/ho)5p'^, we obtain 



doSp = -V&P + i^V^SvL + ( C + ^^7 ) V(V ■ 6vl), (89) 
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which coincides with Eq. (13). We have shown that the hnearized Landau equations, Eqs (77), (83) and (89), are 
derived by choosing Se, Sn, and Sp'' as slow variables. 

Before closing this subsection, let us consider the detail of <?„„(z, k). The memory function can be written as 

<P{z, k) = ~{S{z, k) - iQ{k)S{z, k))- -, (90) 

1 + ^[z, k) 

where we defined 

S{t)^{A^{t),A^), (91) 

S{t) = {iCA^{t),A"^), (92) 

S{t) = {{iCfArXt),A^'). (93) 

For the derivation of Eq. (90), see Appendix B. Since the time derivative of a conserved charge variable is slow, S is 
of order fc; then, we can estimate 1/(1 + S') = 1 + 0{k). And then, <P„„(z, k) becomes 



<Pnn{z,k)^- ^{z,k) - if2{k)^{z,k) 



k ( -jj(z,o)- — ^pj(z,o; 



(94) 



From ^ = — 1 + z^, 



^pj{z,0) = -^gpj = -^noTo, (95) 



where we used ^pj{z, 0) = 0. Therefore, we can write 

^>„„(z,fc) fc2 (^H,,(z,0) - = k^Xiz), (96) 

where we defined the frequency-dependent-thermal conductivity A(z). At z = 0, A(0) coincide with A. This expression 
will be used in the next subsection to derive the linearized Eckart equation. 

B. Linearized Eckart equation 

Next, we derive the linearized Eckart equation. The charge does not dissipate in the Eckart equation, which implies 
that the fluid velocity is chosen as Sj'^/n. Therefore, we choose as a slow variable in addition to {Se, Sp^, Sn}. In 
order to derive the Eckart equation, we first derive the equations of motion for {Se, Sp'', Sn, Sj'}; and then, we remove 
the degrees of freedom of Sp' . 

The equation of time evolution for Se and Sn arc trivial thanks to the conservation laws, 

doSe = -V • Sp, (97) 

doSn = -V ■ Sj. (98) 

As usual, we calculate the reversible terms: 

iQ,p.{k) = J d^xe-'''-''{iCSe{x),Sp'{0)) = -ik'Toho + C'(fc^), (99) 

i%.(fc) = J d^xe-'''-''{iCSn{x),Sf{0)) = -ik'T^n^ + 0{k^), (100) 

if2nf{k) = J cPxe-'^-''{iCSn{x),Sj\<d)) = -ik^gj^^^ + 0{k^). (101) 

Here the explicit form of gjiji is not obtained; however, it is irrelevant in the leading order of the fluid equations, as 
will be seen later. The reversible term for doSp' is the same as that of the Landau equation in Eq. (84). The reversible 
term for doSj' becomes 



i/?j„A" + in^^A" ^ -ik' gjjSiPfj.) - TonoS/3 . (102) 
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Since Sj^ and Sp^ are chosen as slow variables, the inverse metric contain mixing terms: 



gPP gPJ 
gjP gjj 



1 



9j] -9pj \ _ 



9pp933 ~ 9pj9ip \ 93P 9pp 
Then, the conjugate variables for (5j' and 5p^ are 



/3o 



Po9jj -no 
-no /iQ 



/3o 



hoPogjj - 



no 
ho 



A^ 



(103) 

(104) 
(105) 



where we defined the fluid velocity in the Eckart frame as Sv'^^ = Sj^ /no. The memory functions appearing in doSp^ 
are (l>pp and (Ppj, which are both of order fc^. Here, we assume that A^ is of order k, which will be checked later. 
Then, 'PpjA^ can be neglected because it is of order k^. Similarly, <l>ppAP ~ ~/3<PppSvL. Therefore, 



do5p = -VSP + r^V^SvL +[C+\v\ V(V ■ 5vl). 



(106) 



This is the same equation as the Landau one. 

Finally, let us consider the equation for do5j', which is written in the Laplace space: 

znoSv^ = ifijnA"" + ifljeA" - (l>jjA^ - ^jp^^ + no^VEi^ = 0) 

= -g,, VM(/?M) - ^V^-^T - ^,,A^ - <PjpAP + noSv^Eit - 0), 
Jo 



(107) 



where SvE{t = 0) is the initial value of the fluctuation. The important point is that <l>jj is not slow and gives a 
contribution of order k^. This equation will give the relation between these fluid velocities. Since we are interested in 
the first order hydrodynamic equation, it is enough to take into account the difference of the fluid velocity up to of 
order k^. In this order, we can neglect (l>jp because it is of order k^. 
Let us, now, estimate ^jj. From Eq. (90), we obtain 



iS{z,k)~in{k)E{z,k)) 



l + S{z,k) 



-z + 



-zgn 



1 



S{z,k) 
1 



if]{k) 



(108) 



33 



Siz,k) 



33 



where we used S = zS — if2, S ~ zE ~ 1, and if2jj{k) = 0. We may estimate 'Pjj{z,k) at fe = in the leading 
order. First, we consider Sf.pi(z, k). The energy conservation provides zS^pi^z, k) = —iy Spjpi{z, k). Thus, at fc = 0, 
S'epi {z, 0) = 0. Similarly, one can show E^ejiz, 0) = Snp{z, 0) = Snjiz, 0) = 0. Therefore, we may consider the terms 
with Sp"^ or 5f . They are estimated at fc = as 



^pp{z,o) :^pj(z,o; 

^jp{z,Q) :Zjj[z.^Oi 



_ ^ gpp 



9pj 



93P z^jj{z,0) 



1 / hoTo 



noTo 



zUoTo !fil + ^A(z)]' 



(109) 



where we used ^ (z, 0) = —1 + z^ [z, 0) = 0, and Eq. (96). Taking into account the metric, we find 

{giiho - nlTof 1 o , 1 1 
<?,-,(z,0 = V ^ " - —ihgri - nlTo)z = . -z. 

Then, the equation of motion becomes 

nozSvE = -gjjVSiPfi) - ^VST H . {6vl - Sve) - noz{5vL - 5ve) + no5vE{t = 0). 

Jo gi3\{z) 



(110) 



(111) 
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From Eq. (Ill), we obtain 

hoSvL = hoSvE + —g"\{z)(7ioz5vL + gjjV6{/3^i) + ^V(5T - iioSvEit 0) 
no V Jo 



(112) 



= hoSvE - \{z){Toz5vL + VST) + \{z)T^gPP{hQz5vL + V5P) - g"'^^^\{z)5vE{t = 0) 

The third term in the last hne is estimated as h^zSv^ + V(5P = ho6vL{t = 0) + ©(V^) from the equation of motion, 
Eq. (106). Then 

hoSvL = hoSvE - Xiz){TozSvL + VST) + \{z){T^g^^ho5vL{t = 0) - gi^'^^SvEit = 0)) 



^2 (113) 
h^SvE - X{z){Tqz5vl + VST - ToSvLit = 0)) - -^X{z)A\t = 0). 

"0 



This equation is exact for an arbitrary z in the first order of the derivative expansion of spatial coordinate. We note 
that the equations (97), (98), (106) and (113) are those of motion for {Se,Sp^,Sn,6j'^}. These equations are first 
derived in this paper by the projection operator method. We also note that these include the Landau and Eckart 
frames in the view point of the projection operator method. 

Now, let us derive the Eckart equation. We assume that there is no conjugate variable for at the initial time, 
i.e., A^(t = 0) = 0. We also assume that the change of the variable is so slow that z expansion is applicable: 

hoSvL = hoSvE - X{Toz5vl + VST - SvL{t = 0)). (114) 

In the coordinate space, this equation becomes 

hoSvL = hoSvE - X{TodoSvL + VST). (115) 

From this equation, one can confirm = —nog^^Sv]^ — Sv'"^) = n^g^^ {T^d^Sv]^ + ik'^ST)/ho ^ k. Therefore, the 
assumption that A^ is of order k is consistent for slow dynamics. Solving Eq. (115) for hoSv^, we obtain 



hoSvL = — — j^ihaSvE - XVST) 



1 

^ + ^rodo''"''^ (116) 
hoSvE - XiTodoSvE + VST), 



where we dropped doVT in the last line because it is higher order. Inserting Eq. (112) into Eq. (106), we find 

doihoSvE - X{TodoSvE + VST)) = -VSP + j^VHve + + V(V • Sve). (117) 
This equation is equal to the linearized Eckart equation, Eq. (19). 



V. DISCUSSION 



Here, we discuss the results obtained in the previous section. In the derivation of the Eckart equation, we eliminate 
the momentum density from the equations for {Se,Sp^ ,Sn,Sj''}. Then, the Eckart equation has {Se^Sn^Sj^} as the 
dynamic variables, apparently. Nevertheless, the Landau's variables {ScjSp^ ,Sn} describe the slow dynamics, even 
for the Eckart equation. To illustrate this fact, we consider the Onsager reciprocal relation in the Eckart equation. If 
we regard {(5e, Sn, as independent variables, the reciprocal relation is violated. In contrast, if we regard Sp^ as an 
independent variable instead of Sj^ , we will find that the relation is satisfied. This result shows that the slow part of 
(5j' is determined by {Se^Sp^ ,Sn}. 

Also, we discuss the instability of the Eckart equation. The Eckart equation has an unstable mode. However, we 
will find that the equations for {(5e, Jp', Jn, ^j*} do not have the unstable mode and also have the same mode as the 
Landau equation does. We will show that the unstable mode arises from the expansion about 9o in Eq. (116). 
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A. Independent variables of the Eckart equation 

Let us consider the correlations, [dode, 6n) and [dodn, Se) . The important point is that the correlations must satisfy 
the relation 

{do6e,Sn) = {daSn,Se), (118) 

which comes from the time reversal properties of Sn, 6e and the equilibrium state. We note that this relation is 
equivalent to the Onsager reciprocal relation in the linear regime [25-27]. From the Eckart equations, Eqs. (17)-(19), 
the correlations are written as 

{doSn{t,r),Se{t,r')) =-noV ■ {SvE{t,r),5e{t,r')), (119) 

(doSeit, r),Snit, r')) = - /lo V • {Sveit, r),6n{t, r')) 

+ x(v^{ST{t,r),Sn{t,r')) + ToV ■ {dodvEit,r),6n{t,r'))y (120) 

Now, to eliminate the time derivative in the last term of Eq. (120), we use the derivative expansion. Namely, we 
approximate 



AToV • doSvE = ATqV • h^\-V5P+ ...) (121) 
'-V^SP, (122) 







where ... denotes the second and higher order terms about the derivative, like rfV^Sv. Here, we first used Eq. (19) 
for do6vE and next neglect the second-order terms, which finally yields third-order terms. Thus, we obtain 

{doSeit, r), Sn{t, r')) -hoV ■ {SvE^t, r), 6n{t, r')) - A V^S{(3fi){t, r),5n{t, r')) , (123) 

where we also used the Gibbs-Duhem relation, Eq. (85). Let us here estimate the equal-time correlations, ((5t>£;, (5n) , 
{SvE^Se) and [S{(3iJ,),Snj . We assume that the time-reversal property of 6ve (Sj) is odd, according to Eq. (55). 
Then, 6ve is orthogonal to Sn: 

{dvE, Se) = {Sve, Sn) = 0, (124) 

by the time-reversal symmetry. On the other hand, (^S{f3ii),Sn) turns out to be S{r — r') from Eq. (66). Finally, we 
obtain the correlations, 

{doSnit,r),Se{t,r')) =0, (125) 

{doSeit, r), Snit, r')) = V'Sir - r'). (126) 

We see that the Onsager relation, Eq. (118), is violated. 

Next, let us regard Sve as a function of {Se,Sp'',Sn}. Namely, we consider Sp as an independent variables, instead 
of Sj. Then, we now use Eq. (116), which gives the relation between Sj and Sp: 

Sp = hoSvE - X{TdaSvL + VST) 

^h,svE^x(^^ysim, ^'''^ 

where we used the derivative expansion and the Gibbs-Duhem relation in the second line. Solving the above relation 
about Sve, we obtain Sve as the function of {Sn,Se,Sp}: 

1 <c„ , ^^ ^oT^ 



SvE{Sn,Se,Sp) = —Sp + X[ V(5(/3/x). (128) 

ho \ K J 

Here, we note that the time-reversal property of SvE{Se,Sp'^ ,Sn) is not odd because those of Sp and S{l3fi) arc odd 
and even, respectively. Substituting Eq. (128) into Eqs. (119) and (123), we find 

{doSnit, r),Se{t, r')) = V • {Sp{t, r),Se{t, r')) - A (^^^ {SiMit, r), <5e(t, r')), (129) 
{doSe{t, r), Sn{t, r')) ^ {Spit, r), Sn{t, r')) . (130) 
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Here, we note that S{l3fi) is orthogonal to Se by Eq. (65). Then, if we assume the time reversal property of 5p is odd, 
we see that the Onsager relation is satisfied: 



{do6n{t,r),Se{t,r')) = {dnSe{t,r),Sn{t,r')) = 0. 



(131) 



Why cannot we regard 6p as a function of {5n, Se, Sj}7 The reason is that Sj is not essentially slow, and then the 
slow motion of that turn out to be described by the actual slow variables. Thus, the time reversal property differs from 
the original operator at the starting point of Sec. III. Actually, we can show that the projection of Sj on {Sn, Se, Sp} 
gives Eq. (128). The derivation is given in the Appendix C. Namely, the actual expression of Eq. (128) is given as 



(132) 



where V is the projector on {Sn,Se,Sp}, i.e., Sj (Sve) in the Eckart equation is projected one and differs from the 
original operator. In consequence, the slow variables for the Landau frame actually describe the slow dynamics even 
for the Eckart frame. 



B. Modes of the equations for {Se, 5p^ ,Sn, Sj'^} 



Now, we study modes of the equations for {Se,Sp'^,Sn,Sj^}, Eqs. (97), (98), (106) and (113). In the Fourier space, 
the equations are written as the following matrix form: 



M 



where we decomposed Sp and Sj into the longitudinal and transverse components: 



/ Se{^,k) \ 






Sn{uj, k) 







Sp\\{uj,k) 














Sp± (w,fc) 







\SjA^,k)J 







Also, we introduced the matrix 

/ - 



M 











k ■ Sp, 
k ■ Sj, 



ik 




Sp± = Sp — Sp^^k, 
Sj± = Sj - (5j||fc. 



ikape ikapn—iuj + k^F^ 
ikikaTeikaxn 



-lUJ - 





ik 


'{ho/nQ)rj 











-iu! + fc2(7]//lo) 










-{ho/no)rj) 



where 



an 



OiTe 



dP 



1 



dP 

dn 



dT 
dn 



(133) 



(134) 
(135) 



(136) 



We here neglected the frequency dependence of the thermal conductivity. We can obtain dispersion 
det M = 0. Those are given as the following to second order in k: 



UJ ' 



'ik^Pt, 
-ik'^{i]/ho). 



(137) 

(138) 

(139) 
relations from 

(140) 
(141) 

(142) 
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where we introduced the thermal and sound diffusion constants, and the sound velocity: 

A 



tlQCp 



dP 

9e J ,/„ 



(143) 
(144) 
(145) 



Here, cp is the specific heat at constant pressure. From Eqs. (140)-(142), we see that the equations for {(5e, (5p*, (5n, 

do not have unstable modes. Moreover, these dispersions are the same as for the Landau equation [34], i.e., the current 

density bf' is irrelevant for slow dynamics. 

When does the unstable mode arise in the derivation of the Eckart equation? The unstable mode arises from the 
expansion about 9o in Eq. (116). Let us now study the modes after the expansion. For simplicity, we consider only 
the transverse components. Those are given as 



-iuj + k'^{ri/ha) 
no/ ho 







■lOJ ■ 



and dispersions are 



UJ = ~ik'^{r]/ho), 



(146) 



(147) 
(148) 



We see that there are two modes; the first mode is the same as Eq. (142) and the second one is unstable. We note 
that, before the expansion, the transverse components have the only one mode. Namely, by the do expansion, the 
number of modes increase and the unstable mode arises. 

Now, we note that the eigenvalue of the unstable mode is not proportional to the wavenumber k and large. In other 
words, the unstable mode is massive and not a hydrodynamic mode. Therefore, we can get rid of the unstable mode 
by introducing a cutoff in the frequency space because hydrodynamics is an effective theory in the low-frequency and 
-wavenumber region. 



VI. SUMMARY 



We studied relativistic hydrodynamics in the linear regime with Mori's projection operator method. From the 
view point of the projection operator method, we discussed that the difference of the frames is not the choices of the 
reference frames but rather those of the slow variables. We also found that the the slow variables for the Landau 
frame are the conserved charges whereas those for the Eckart frame include the current of the conserved charge, 
which is not essentially slow. In fact, we derived the slow dynamics by the projection operator method for the sets, 
{5e, Sp'^ , 6n} and {6e,6p'',5n,6j'^}, as the slow variables. We first showed that the natural choice, Eq. (23), gives the 
linearized Landau equations. Next, we derived the equations of motion for {Se, Sp^, 5n, which include the Landau 
and Eckart frames in the view point of the projection operator method. And then, we derived the linearized Eckart 
equation by eliminating 5p^ from the above equations. 

We also discussed properties of the derived equations. In particular, by considering the Onsager relation, we 
illustrated that the slow part of the particle current, which is the fluid velocity for the Eckart frame, is determined 
by Landau's variables, {5e, 6p^, Sn}. Moreover, we studied how the unstable mode of the Eckart equation arises. We 
found that the do expansion for the elimination of Sp^ causes the unstable mode. Then, we pointed out that the 
unstable mode is "massive" and is ruled out by the cutoff in the frequency space. Also, we found that the equations 
for {Se,Sp^,Sn,Sj^} has the same modes as the Landau equation has. Thus, the particle current Sj' is irrelevant 
for the first order hydrodynamics. Recently, some atuhors also point out that the Landua frame is natural for the 
relativistic hydrodynamics, based on the renormalization group method [32]. 

Here, we note that this study is first for the derivation of relativistic hydrodynamics based on the projection operator 
method. We stress that our derivation is independent of the microscopic details; we determine the metric from the 
Lorentz symmetry on microscopic scale and the thermodynamics. The earlier studies [7-10] assume the relativistic 
Boltzmann equation as the underlying microscopic theory, which is, however, only valid for a weakly-correlated system 
like a dilute gas. Furthermore, we note that our study is independent of what are the local equilibrium and the local 
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rest in the relativistic fluids because we considered the fluctuations from globahy equihbrium state. Instead, our study 
is restricted in the hnear regime. 

Now, we comment on the Lorentz covariance of hnearized hydrodynamics. The hnearized Landau equations, (11)- 
(13), are not Lorentz covariant. The reason is the following: We here considered the fluctuations in the background 
medium. For such a system, the Lorentz transformation boosts the fluctuations but does not the background. Then, 
the boosted system differs from that before the boost. Thus, the linearized equations are valid only for the rest 
frame of the medium and not Lorentz covariant. Actually, by the same reason, the Navier-Stokes equation is Galilei 
covariant whereas the linearized one is not covariant. However, we note that the linearized Navier-Stokes equation 
well describes the fluctuation in a fluid [17, 28]. Namely, linearized hydrodynamics has no problem in the linear 
regime. 

In this paper, we used the linear projection operator, and our study is restricted in the linear regime. Then, 
it is interesting to derive relativistic hydrodynamics by the nonlinear projection operator [29]. We note that, for 
nonrelativistic fluids, the full Navier-Stokes equation is derived by the nonlinear one [30]. Furthermore, we here 
focused and discussed the Landau and Eckart frames. Then, we concluded that the Landau frame is natural for the 
slow dynamics. However, other open problems of relativistic hydrodynamics, such as the acausal propagation, have 
not been solved yet. It would be interesting to discuss these in the view point of the projection operator method. The 
projection operator method may give an insight into these problems, as well as that of the frames. 
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Appendix A: Properties of inner product 

Here, we summarize properties of the inner product deflned in Eq. (28). For Hermitian operators, 

{A{t),B{0)) = {A{t),BiO)r = (5(0), A{t)), (Al) 
iiCA{t), B{0)) = -(Ait), iCBiO)), (A2) 

{tCA{t), B{0)) = -^{[A{t), i?(0)])eq, (A3) 

{A{t),B{0)) = eAeB{A{-t),B{0)) = eAeB{B{t), A{0)) (A4) 

are satisfled. Here sa a-nd eg denote the sign associated with time reversal transformation, which is deflned with the 
time reveal operator T as T^^A{t)T = eAA{—t). 



Appendix B: Memory function 



In this Appendix, we derive the full expression of the memory function following Ref. [31]. Let us flrst deflne the 
key functions: 

S{t) = {tCA,,{t),A^), (Bl) 

^(t)^((z/:)2A„(t),A™). (B2) 

The reversible term can be expressed by if2 = S{0). S{t) coincides with the memory function obtained by replacing 
Q with unity. Expanding exp(tQi£) in terms of ViC, we obtain 

e*2'^ = e**^(l + V(-l)" / dil--- / dt„z/:^(ti)---i£^(t„)), (B3) 

where 

C^{t) EE e-'^'VCe'^*. (B4) 
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Acting C^{t) to An; wc obtain 

iCP{t)Ara{-t') = e-'^'PiCe'^'A,n{-t') = An{^t){iCA,n{t - t'), A") = - t')Ani-t). 

Then, for an operator O, 



is satisfied. Here we defined 



(l + V(-l)"/ dh--- dt„lCPih)---lC''itn))0 
^ n=l ^ 

X dp{tn)S(tn-l - tn)S{tn-2 - tn-l) ■ ■ ■ ^{tl - t2)A{t - h) 

dp{t) = {iCO{t),A"'). 



Performing the Laplace transform, we obtain 



1 



(B5) 



(B6) 



(B7) 



(B8) 



Using this equation, we find tlie full expressions of the nose term and the memory function in the Fourier-Laplace 



space: 



k) = iCA{z, k) - {iQ + E:{z, k)) 



<Z>(z, fc) = -{E:{z, k) - iQ{k)E:{z, k)) 



1 



1 + S(z,fc) 
1 

l + ^(z,fc)' 



A(z,fc), 



(B9) 
(BIO) 



Appendix C: Projection of 5j on {5e,5p,5n\ 
Here, we show that the projection of 5j on {5n, Se, dp} gives Eq. (128). Namely, we will show 



P 



Sj{t,x) - !I^Jp{t,x) + x(^ ] V6{l3fi){t,x)] 



to the first order in k. In the Fourier-Laplace space, Eq. (CI) becomes 

P \sj{z, fc) - ^5p(z, fc) + ^k\('-^) 5(/3/.)(z, fc)] 
ho V "0 / 

Here we decompose Eq. (C2) into the longitudinal and transverse components: 



P 



Sj\\{z,k) - ^5p||(z,fc) + ifcA^^^^^ S{Pn){z,k)] 



P 



Sjj_{z,k) - '!^5p±{z,k) 
no 



= 0, 

= 0. 

= 0, 
= 0. 



(CI) 

(C2) 

(C3) 
(C4) 



Let us show the equation for the transverse component, Eq. (C4). Consider the projections, P5j^{z;k) and 
PSp±{z, fc), which are given as 



PSj±{z, fc) = S'j-^"(z, fc)(5n(0, fc) + Sj^^iz, k)Se{0, fc) -t- Sj^p^{z, k)Sp±{0, fc), 
PSp_l{z, fc) = Sp^^iz, k)Sn{0, fc) + Ep/{z, k)5e{Q, fc) + Sp/^{z, k)Sp±{0, fc). 



(C5) 
(C6) 



We note that fc expansions of Sj^"{z,k) and Ejj^'^{z,k) give only odd-order terms in fc from the tensor structure. 
Then, we can drop Sj^^'{z,k) and Sj^'^{z,k) because the odd terms are orthogonal to the transverse component. 
Then, Eq (C4) turns out to be 



(5px(0,fc) = 0. 



(C7) 
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Here, let us consider Sj^^-^{z,k) and Sp^'^^ {z,k). For this task, we now use the equations of motion for 
{6n,6e,6p,5j}, Eqs. (97), (98), (106) and (113). From these equations, we obtain the equations for the transverse 
components as 



z + Pri_ 

z + Fj -{ho/no)rj 



5p±{z,k)\ _ 
Sj±{z,k) 



Spj_{t = 0,k) 

6p^[t ^Q,k) + noTog^^ (sj^it = 0, fc) - 7i^'5pj.(t = 0, fc) 



Therefore, 



dp±iz,k) 
dj±{z,k) 



1 



z + k'^ri_ 

{no/ho) 



dp^{t^O,k), 



5p^{t = 0,A;) 



z + fc2rj_ 

If we notice that the conjugate variable for Sp^ is 
in the space of {5e, Sp, Sn}, we find 



hoFj 



hoTo 



■5p\ 



^Vi^^ (z, fe) — 



1 



z + k'^r^ ' 



where we used Eqs. (67) and (68). Finally, we arrive at 



tlQ 



no 



5j^{t = 0, fe) - T^5p±(t = 0, fe 



(C8) 

(C9) 
(CIO) 

(Cll) 

(C12) 
(C13) 

(C14) 



We have shown the equation for the longitudinal component, Eq. (C4). By the similar procedure, we can show that 
for the longitudinal component, Eq. (C3). 
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